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Abstract 

Prompted by recent solar and atmospheric data, we re-analyze the four-neutrino 
Q_i' description of current global neutrino oscillation data, including the LSND evidence 



for oscillations. The higher degree of rejection for non-active solar and atmospheric 
oscillation solutions implied by the SNO neutral current result as well as by the latest 
' 1489-day Super-K atmospheric neutrino data allows us to rule out (2-1-2) oscillation 

schemes proposed to reconcile LSND with the rest of current neutrino oscillation 
^ ! data. Using an improved goodness of fit (g.o.f.) method especially sensitive to the 

?H I combination of data sets we obtain a g.o.f. of only 1.6 x 10~^ for (2-1-2) schemes. 

Further, we re-evaluate the status of (3-1-1) oscillations using two different analyses 
of the LSND data sample. We find that also (3-1-1) schemes are strongly disfavoured 
by the data. Depending on the LSND analysis we obtain a g.o.f. of 5.6 x 10^^ or 
7.6 X 10~^. This leads to the conclusion that all four-neutrino descriptions of the 
LSND anomaly, both in (2-1-2) as well as (3-1-1) realizations, are highly disfavoured. 
Our analysis brings the LSND hint to a more puzzling status. 
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Figure 1. The six types of four-neutrino mass spectra. The spacings in the vertical 
axes correspond to the different scales of mass-squared differences required in solar, 
atmospheric and short baseline neutrino oscillations. 

1 Introduction 



The atmospheric neutrino data |l|, ||, ||, including the most recent 1489 Super- 
K data sample provide strong evidence for u^j, oscillations into an active neu- 
trino (mainly z/^) with Am^^^ ~ 2 x 10~^ eV^ ||^. On the other hand, apart 
from confirming, once again, the long-standing solar neutrino problem |^, 
1^, H, 0, the recent results from the Sudbury Neutrino Observatory (SNO) 
|1T0| have given strong evidence that solar neutrinos convert mainly to an ac- 
tive neutrino flavor. This suggests that an extension of the Standard Model 
of particle physics is necessary in the lepton sector, capable of incorporating 
the required 



11, 1^ 



the 



conversion. Although certainly not yet unique 
most popular explanation of the solar neutrino data is provided by the active 
large mixing angle (LMA) neutrino oscillation hypothesis, characterized by a 
neutrino mass-squared difference ArrigQ^ ^ 10""^ eV^ M. 



In contrast, reactor and accelerator neutrino data |T^, 0, |T^, [16], [I^ lead to 



negative results. However, the LSND experiment [|T^, has provided positive 
results, which may or may not be confirmed by the forth-coming MiniBooNE 

> 



20 1 . The required neutrino mass-squared difference Ami 



experiment 

0.2 eV^ is in conflict with the neutrino mass-squared differences indicated by 
solar and atmospheric data in a minimal three-neutrino picture. Four-neutrino 
models |22|, |2^, ^ ^ ^ potentially account for all current oscillation 
data. The status of four-neutrino descriptions has been presented in Ref 
An exhaustive list of four- neutrino references can be found in Ref. |29]. 
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The six possible four- neutrino mass spectra are shown in Fig. 1. For the case 
Am^gt^n ^ ^"'-ATM' ^"^soL' which we tacitly assume in order to reconcile the 
LSND anomaly with solar and atmospheric data, these schemes can naturally 
be divided into two very different classes, usually called (3-1-1) and (2-1-2) mass 
schemes [^]. It is important to note that (3-1-1) mass spectra include the three- 
active neutrino scenario as limiting case. In this case solar and atmospheric 



2 



neutrino oscillations are explained by active neutrino oscillations, with mass- 
squared differences ArngQ^ and Am^^^j, and the fourth neutrino state gets 
completely decoupled. We will refer to such limiting scenario as (3+0). In 
contrast, the (2+2) spectrum is intrinsically different, as there must be a 
significant contribution of the sterile neutrino either in solar or in atmospheric 
neutrino oscillations or in both. 



In recent studies |^ it has been realized that there is considerably tension 
in four-neutrino fits to the global data. In such global four-neutrino analyses 
one is faced with the problem that there are different data sets, which all give 
very good fits if analyzed separately. Problems arise due to the combination 
of the different data sets in a four-neutrino framework. In this letter we re- 
evaluate the status of four-neutrino interpretations of the LSND anomaly in 
the light of the recent solar [|, |^, as well as atmospheric [0 neutrino 
results. To evaluate the quality of the fit we will apply appropriate statistical 
methods, which are especially sensitive to the combination of different data 
sets in a global analysis. We find that both the SNO NC result as well as the 
1489-day Super-K atmospheric neutrino data strongly reject against sterile 
neutrino conversions. This essentially rules out (2+2) mass schemes. We also 
re-evaluate the status of (3+1) schemes by considering two different analyses 
of LSND data [1^, pl|]. We will further elaborate the result of previous studies 



30| , |33| , [3^ , p5| , |36| , that in (3+1) schemes the LSND signal is in serious 



disagreement with bounds from short-baseline (SBL) experiments reporting 
no evidence for oscillations IIT3|, [T^, Ol . The net result is that all four-neutrino 



descriptions of the LSND anomaly, both in (2+2) as well as (3+1) realizations, 
are strongly disfavoured by the data. This brings the LSND anomaly to a 
theoretically more puzzling status. We note that also cosmology put strong 
constraints on four- neutrino schemes (for recent analyses see Ref. [B3]). 



The plan of the paper is as follows. In Sec. 2 we briefly describe our param- 
eterization and approximations used in the global four-neutrino analysis. In 
Sees. 3 and 4 we summarize the solar and atmospheric neutrino data and their 
analysis Q. In Sec. 5 we describe the SBL data we are using, and we com- 



pare the two different analyses of the LSND data |T^, . In Sec. 6 we show 
that (2+2) oscillation schemes are ruled out because of the tension between 
solar and atmospheric data, whereas in Sec. 7 we present our re-analysis of 
the disagreement in SBL data in (3+1) oscillation schemes in light of the two 
LSND samples. In Sec. 8 we show that the goodness of fit is very bad in all 
four-neutrino cases. Furthermore, we compare the relative quality of the fit for 
the schemes (3+1) and (2+2), as well as the three active neutrino case (3+0). 
The quantitative statistical criteria to analyze the conflict between different 
data sets, which are used in Sees. 6, 7 and 8, are formulated in the appendix A. 
In summary, we find that all four-neutrino descriptions of the LSND anomaly 
both in (2+2) as well as (3+1) realizations are highly disfavoured due to recent 
data, as mentioned in our conclusions, Sec. 9. 
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2 Four-neutrino oscillations 



In this letter we interpret current neutrino oscillation data in a four-neutrino 
framework. This requires the existence of a fourth light neutrino, which must 



be sterile to account for the LEP invisible Z width data |2^, ^ |29[| . 

The simplest theoretical models add just one such hght SU{2) (g) U{1) singlet 
neutrino and ascribe its lightness either to an underlying protecting symmetry, 
such as lepton number |22[, or to the existence of extra dimensions p5| , |26[| . 
Symmetry breaking is required in order to generate solar and atmospheric 
oscillations, characterized by splittings which arise due to small tree-level 
p7| or radiative effects [^, Many other four- neutrino models have been 



considered ||2^, . 



Quite generally, the four-neutrino charged current leptonic weak interaction 
is specified as by a rectangular 3x4 lepton mixing matrix K = QU, with Q 
diagonalizing the 3x3 charged lepton mass matrix and U diagonalizing the 
4x4 Majorana neutrino mass matrix. Altogether the elements of the mixing 
matrix K are characterized by 6 mixing angles and 3 physical phases ||39|| , 
which could lead to CP violation in oscillation phenomena For simplic- 



ity we neglect CP-violating phases, whose effects are expected to be small in 
the experiments we consider, due to the strong hierarchy of the mass-squared 
differences required by the experimental data. This leaves us with nine param- 
eters altogether relevant for the description of CP conserving four-neutrino os- 
cillations: 6 mixing angles contained in the mixing matrix and 3 mass-squared 
differences. 



As shown in Ref. it is very important to adopt a convenient choice for these 



parameters in order to factorize the analysis. Here we adopt all the notation 
and conventions introduced in Ref. PS[, which are based on the use of physi- 



cally relevant quantities. The 6 parameters ArngQ^; ^sol, ^t^atm^ ^atm, ^"^lsnd' 
^LSND are similar to the 2-neutrino mass-squared differences and mixing angles 
and are directly related to the oscillations in solar, atmospheric and the LSND 
experiments. For the remaining 3 parameters we use rjs, rje and d^. Here, rjs {rje) 
is the fraction of fg {^e) participating in solar oscillations, and 1 — rf^ is the 
fraction of participating in atmospheric oscillations (for exact definitions 
see Ref. [0). For the analysis we adopt the following approximations: 

• We make use of the hierarchy AmgQL ^ ^"^atm ^ ^'^Lnd- This means that 
for each data set we consider only one mass-squared difference, the other 
two are set either to zero or to infinity. 

• In the analyses of solar and atmospheric data (but not for SBL data) we 
set r^e = 1, which is justified because of strong constraints from reactor 
experiments M, 0, ITBI. 
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Figure 2. Parameter dependence of the different data sets in our parameterization. 

Due to these approximations (well justified by reactor data and in oder to 
fit in LSND with solar + atmospheric data) the parameter structure of the 
problem gets rather simple The parameter dependence of the different 
data sets solar, atmospheric and SBL is illustrated in Fig. 2. We see that only 
rjs links solar and atmospheric data and links atmospheric and SBL data. 
All the other parameters are "private" to one data set. 



3 Solar neutrinos 



Here we borrow from Ref. 0] a recent re-analysis of present solar data in- 
cluding a study of degree of rejection against the sterile neutrino oscillation 
hypothesis. The experimental data used are the solar neutrino rate of the chlo- 
rine experiment Homestake 0, the rates of the gallium experiments SAGE 0, 
GALLEX and GNO , as well as the 1496-day Super-Kamiokande data sam- 
ple 1^ in form of 44 recoil electron zenith-energy bins. Moreover, we include 
the data from the recent SNO charged current and neutral current event mea- 
surements, as well as day-night information in the form of 34 data points of 
the observed energy spectrum For the solar neutrino fluxes we use the 



Standard Solar Model fluxes M], including its prediction of the B flux. The 



first hints in favour of the LMA solution |^ , which followed mainly from the 
flatness of the Super-K spectrum, have now become a robust result, thanks to 
the new SNO data and additional Super-K spectrum data. In contrast with 
the pre-SNO-NC situation, presently all non-LMA solutions are strongly re- 
jected, so that in this paper we will restrict the solar oscillation parameters 
to the LMA region. 

Another important consequence of the recent SNO-NC data is specially rel- 
evant for us, namely, that the degree of rejection against solar neutrino os- 
cillations to sterile neutrinos has now become much more significant. This 
results can be summarized in Fig. 3. The line labeled "after" shows ^Xsol 
a function of < 77s < 1, after recent solar neutrino up-dates, most important 



SNO-NC |]T0[ . We minimize with respect to the other two solar neutrino oscil- 
lation parameters Am^Q^ and ^sol, required only to lie in the LMA region. For 
comparison we show also the line labeled "before" , corresponding to the data 
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Figure 3. A^sol foi' the favoured LMA solution, before and after SNO NC data |To[| , 
displayed as a function of r^^- 



used in Ref. p8[. From this figure one can place a mass scheme- independent 
limit on the parameter r\s from solar data: 



solar data (LMA): 77^ < 0.45 



C.L.I 



compared to ris < 0.52 obtained in Ref. p3]. Moreover, the pure sterile solution 
{ris = 1) and large values of t]s are now much more disfavoured than previously. 



4 Atmospheric neutrinos 



For the atmospheric data analysis |^ we use the following data from the Super- 
Kamiokande experiment [0]: e-like and /i-like samples of sub- and multi-GeV 
data, as well as the up-going stopping and through-going muon data. The 
Super-K data correspond to the 1489-day sample presented at Neutrino 2002 
0. In contrast to previous analyses p8| , ^ ^ here we adopt a more precise 
ten-bin presentation of the corresponding contained zenith-angle distributions. 
Further, we use the data from MACRO |]^, including the recent update of their 
up-going muon sample (10 angular bins). Since ArngQ^ ^"^atm adopt 
the standard approximation of neglecting Aml^^ when performing the analysis 
of the atmospheric dataQ. For further details of the atmospheric neutrino 
analysis see Refs. E3, H5| and references therein. 



Let us consider the impact of atmospheric data alone on the parameter rjg 



As shown in Ref. [2S] the behaviour of the atmospheric data as a function of 



For detailed studies of this issue see Ref. 145 
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Figure 4. Ax^tm displayed as a function of rjs for (3+1) and (2+2) schemes, before 
and after the most recent Super-K up-date Q. 

ris depends on the mass schemes considered, with three inequivalent casesQ 
(3+l)a,d, {S+l)b,c and (2+2)^,5- The quaUtative behaviour of (3+l)a,d and 
(3+l)b^c schemes (see Fig. 1) is similar, however, small differences appear due 
to matter effects in atmospheric oscillations. In the following we will always 
minimize between these two cases and refer simply to the (3+1) scheme. 



We present in Fig. 4 the ^Xatm profile displayed as a function oi < rjs < 1 
minimized with respect to all other parameters, before and after the latest 
data of Ref. have been included. We see that the new data significantly im- 
prove the sensitivity of atmospheric data with respect to a sterile component. 
Especially for (2+2) schemes small values of 7]^ are highly disfavoured, which 
will be crucial for the exclusion of (2+2) spectra. From the figure we obtain 
the following bounds on r]s at 99% C.L.: 



atmospheric data: 



rjs < 0.32 for (3+1) schemes, 
Tjs > 0.65 for (2+2) schemes. 



5 Reactor and accelerator data 



We divide the data from reactor and accelerator experiments in data from 
experiments reporting no evidence (NEV) for oscillations and in the LSND 
experiment. In this paper we compare the implications of two different analyses 



^ Note that rjg is defined as the fraction of sterile neutrinos in solar oscillations 
and is only indirectly related to the fraction of sterile neutrinos in atmospheric 



oscillations. This relation is different for the various mass schemes |28|. 
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Figure 5. Comparison of the 99% C.L. regions of the global LSND analysis ll| 
(shaded region) and of the analysis performed in Ref. ||3l| (solid line) using mainly 
the decay-at-rest (DAR) data sample. 

of the LSND data. First, we use the likelihood function obtained in the final 



LSND analysis [19] from their global data with an energy range of 20 < 



Ep < 200 MeV and no constraint on the likelihood ratio R-y (see Ref. IT 



for details). This sample contains 5697 events including decay-at-rest (DAR) 
i^/^ i^e, and decay-in-flight (DIF) z/^ ly^ data. We refer to this analysis 
as LSND global. Second, we use the LSND analysis performed in Ref. |5T[ 
based on 1032 events obtained from the energy range 20 < E^, < 60 MeV and 
applying a cut of > 10~^. These cuts eliminate most of the DIF events 
from the sample, leaving mainly the DAR data, which are more sensitive to 
the oscillation signal. We refer to this analysis as LSND DAR. 



In both cases the collaboration has provided us with the likelihood function 
obtained in their analyses [|^, ^ which we convert into a according to oc 
— 21n£ (see Ref. |^ for details). In Fig. 5 we compare the 99% C.L. regions 
obtained from the two LSND analyses. The LSND DAR data prefers somewhat 
larger mixing angles, which will lead to a stronger disagreement of the data in 
(3+1) oscillation schemes. Furthermore, the differences in between the best 
fit point and no oscillations for the two analyses are given by Ax^^^g^, = 29 
(global) and A^^^^g^ = 47 (DAR). This shows that the information leading to 
the positive oscillation signal seems to be more condensed in the DAR data. 

Concerning the NEV data, we take into account the detailed information from 
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the SBL disappearance experiments Bugey and CDHS [0, as described 
in Ref. |3^. The constraints from KARMEN [13| are included by means of 
the hkehhood function provided by the collaboration Further, we include 
the results on the survival probability from the long-baseline experiments 
CHOOZ [g (P = 1.01 ± 0.028 ± 0.027) and Palo Verde 
0.024 ± 0.053) as explained in Ref. l2^ 



m {P = 1.01 ± 



6 (2+2) oscillations: ruled out by solar and atmospheric data 



The strong preference of oscillations into active neutrinos in solar and at- 
mospheric oscillations leads to a direct conflict in (2+2) oscillation schemes 
p8| , |3^ , This is evident from the comparison of Eqs. (1) and (2). We will 
now show that thanks to the new SNO solar neutrino data |TD[ and the im- 
proved SK statistic on atmospheric neutrinos 0] the tension in the data has 
become so strong that (2+2) oscillation schemes are, in contrast to Ref. ||28[| , 
essentially ruled out. 

In Fig. 6 we show the a function of r]s for solar data and for atmospheric 

combined with SBL data[3]. Furthermore, we show the of the global data 
defined by 

x\Vs) = ^xloM + AXatm+sbl(^s) • (3) 

In this section we investigate quantitatively the degree of disagreement of the 
two data sets solar and atmospheric+SBL in (2+2) oscillation schemes using 
the statistical methods described in detail in the appendix. 

First, in the method which we are calling parameter consistency test (PC) we 
calculate the CL. at which still some value of the common parameter rjs is 
allowed by both data sets. From Fig. 6 we find 

Xlc = 12.2. (4) 

This implies that only if we take both data sets at the 99.95% CL. a value of 
rjs exists, which is contained in the allowed regions of both sets. 

Second, the parameter goodness of fit (PG) makes use of the defined in 
Eq. (3). This criterion evaluates the g.o.f. of the combination of data sets, 
without being diluted by a large number of data points, as it happens for 



^ For completeness we use here (atmospheric+NEV+LSND global) data. However, 
the impact of SBL data on the x^-dependence on rjs is very small, as can be seen 
by comparing the corresponding lines in Figs. 6 and 4. Also the differences between 
LSND global and LSND DAR are negligible in this case. 
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Figure 6. AxioD ^xItm+sbl and Xgiobai as a function of % in (2+2) oscillation 
schemes. Also shown are the values Xpc and Xpc relevant for parameter consistency 
and parameter g.o.f., respectively. 

the usual g.o.f. criterion. For the case of two data sets constraining one single 
common parameter this method is equivalent to subtracting the best fit values 
for the parameter obtained from the two data sets and calculating the CL. at 
which the difference is consistent with zero. According to Eq. (A. 5) we find 



-2 



X 



23.5. 



(5) 



Evaluated for 1 d.o.f.Q such a high x^-value leads to the marginal parameter 
g.o.f. of 1.3 X 10"^. 

From Eqs. (4) and (5) we conclude that (2+2) oscillation schemes are highly 
disfavoured, if not ruled out: solar data and atmospheric+SBL data are in 
disagreement at the 3.5a level according to the PC test and at the 4.8cr level 
according to the PG. This is a very robust result, independent of whether 
LSND is confirmed or disproved. We note, however, that in the unlikely event 



that the KamLAND reactor experiment |^ will see no signal and rule out the 
LMA solution, the fit of (2+2) schemes may improve. In that case, although the 
solar fit quality will substantially deteriorate since we are left with non-LMA 
solar solutions, the agreement with atmospheric data will improve, especially 
for the quasi- vacuum case M. 



Applying Eq. (A. 5) we find for the number of d.o.f. SgoL + 7atm+sbl 
corresponding to the single common parameter r]s. 
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7 (3+1) oscillations: strongly disfavoured by SBL data 



It is known for a long time that (3+1) mass schemes are disfavoured by the 



comparison of SBL disappearance data with the LSND result [p2| , p3| . Here we 
extend the first quantitative analyses of Refs. |^ ^ by evaluating whether 
the two LSND data sets (global [|l^] and DAR ||3T|) are consistent with all the 
rest of the global sample of oscillation dataQ (SOL+ATM+NEV) by using 
the statistical methods explained in detail in the appendix. 



Following previous works |3^, |3^ and the original analyses of Bugey [|1J] and 



CDHS [0 we perform the analysis for fixed values of Aml^j^-^. Therefore, the 
two data sets we are comparing depend on Ilsnd and 8sol+atm+nev param- 
eters and have only the single parameter ^lsnd in common. For each value of 
Am^gt^n the analysis is completely analogous to the one of the previous sec- 
tion. In Fig. 7 we show the results of the PC test and the PG test for both 
LSND data sets as a function of Amlg^^. The horizontal lines indicate the 
C.L. of 95.45%, 99.73% and 99.9937%, corresponding to 2, 3 and 4 standard 
deviations, respectively. From this figure we observe that the disagreement is 
roughly one standard deviation worse with the LSND DAR analysis compared 
to the LSND global data analysis. The larger values of sin^ 26'lsnd preferred 
by the DAR data are in stronger disagreement with the bounds of the NEV 
experiments. 

The upper panel of Fig. 7 shows that for the global LSND data a value of 
^LSND consistent with the rest of the oscillation data exists at best at the 2a 
level. For LSND DAR parameter consistency is possible only at a C.L. greater 
than 99% (for most Am^gND values, only at more than 3a). The parameter 
g.o.f. is shown in the lower panel. For the global LSND data it is always worse 
than 3a, except for Amlgj^^ ~ 6 eV^ it reaches 1%. For LSND DAR it is for 
most values of Aml^j^^ below 0.01%. The weakest disagreement in all cases 
occurs for the "LSND island" around Aml^j^j^ ~ 6 eV^. Note that this does 
not correspond to the best fit value in a global fit (see next section). For the 
LSND global data the best fit point is at Aml^^^ = 27 eV^ with a PC at 97% 
C.L. and a PG 0.2%. For LSND DAR the best fit point is at Aml^^^ = 0.91 
eV^ with a PC at 3a and an PG of 3 x 10"^ 



^ Here the NEV data plays a key role, but also atmospheric data is important in 
the lower regions of Am^gNQ |33, 37 1. For completeness we include also solar data, 
although the impact on the LSND amplitude is very small. 
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Figure 7. Compatibility of LSND with solar+atmospheric+NEV data in (3+1) 
schemes. In the upper panel we show the C.L. of the parameter consistency whereas 
in the lower panel we show the parameter g.o.f. for fixed values of AttIlsnd- The 
analysis is performed both for the global [19| and for the DAR |31| LSND data 
samples. 

8 Comparing (3+1), (2+2) and (3+0) hypotheses 



With the methods developed in Ref. ||2^ we are able to perform a global fit to 



the oscillation data in the four-neutrino framework. This approach allows to 
statistically compare the different hypotheses. Let us first evaluate the g.o.f. 
of (3+1) and (2+2) spectra with the help of the PG method described in the 
appendix. We divide the global oscillation data in the 4 data sets SOL, ATM, 
NEV and LSND. According to Eq. (A. 5) we consider the x^-function 



+ Axj,gy(^LSND) AmLSND) dfj,, rjg) + Ax^sndI^lsndj AmLSNp) , 



(6) 
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Table 1 

XpQ for (3+1) and (2+2) oscillation schemes in the global analysis for the two LSND 
data sets global |jl^ and DAR [^]. In the last row the parameter g.o.f. is shown for 
4 d.o.f.. 

and the relevant number of d.o.f. is given by 



3sOL + 4aTM + 4nEV + 2lSND — Qparameters — 4 . (7) 

This corresponds to the 4 parameters ris,d^, Olsnd, ^"^lsnd describing the cou- 
pling of the different data sets. 



The results of this analysis is shown in Tab. 1. We show the contributions of 
the 4 data sets to Xpg = xLin for (3+1) and (2+2) oscillation schemes and 
for both LSND analyses. As expected we observe that in (3+1) schemes the 
main contribution comes from SBL data due to the tension between LSND 
and NEV data in these schemes. For (2+2) oscillation schemes a large part 
of XpG comes from solar and atmospheric data, however, also SBL data con- 
tributes significantly. This comes mainly from the tension between LSND and 
KARMEN which does not depend on the mass scheme and, hence, also 
contributes in the case of (2+2). Therefore, the values of Xpg given in Tab. 1 
for (2+2) schemes are higher than the one of Eq. (5), where the tension in 
SBL data is not included. 



In the last row of Tab. 1 we show the PG for the various cases. The best g.o.f. 
is obtained with the global LSND data in the (3+1) scheme. However, even in 
this best case the g.o.f. is only 0.56%. It is fair to say that all other cases are 
essentially ruled out by the PG criterion. In general the fit is worse with the 
LSND DAR analysis. Due to the larger values of sin^ "iOi^sm preferred by this 
data the disagreement between LSND and NEV experiments is significantly 
stronger than with the global data. We observe that the fit gets worse, if one 
selects LSND data with higher sensitivity to the oscillation signal. This seems 
rather discouraging for the four-neutrino hypothesis as an interpretation of 
LSND. 
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Ax' 


(2+2) (3+0) 


LSND global 
LSND DAR 


17.8 20.0 
11.5 28.5 



Table 2 

Ax^ of (2+2) and (3+0) oscillation schemes with respect to (3+1) spectra for the 
two LSND data sets global [|l9| and DAR [g^]. 

Although we have seen that none of the four-neutrino mass schemes can pro- 
vide a reasonable good fit to the global oscillation data including LSND, it 
might be interesting to consider the relative status of the three hypotheses 
(3+1), (2+2) and the three-active neutrino scenario (3+0). This can be done 
by comparing the x' values of the best fit point (which is in the (3+1) scheme) 
with the one corresponding to (2+2) and (3+0). In Tab. 2 we show the Ax' 
of (2+2) and (3+0) with respect to (3+1) for LSND global and DAR data. 

First we observe that for LSND global (2+2) is strongly disfavoured with 
respect to (3+1) with a Ax' = 17.8. The reason for the big change with respect 
to the value of Ax' = 3.7 found in Ref. [^] is the improved sensitivity of 
solar (SNO NC) and atmospheric (SK 1489-days) data to a sterile component. 
With this new data now (3+1) schemes are clearly preferred over (2+2): for 
4 d.o.f. a Ax' = 17.8 implies that (2+2) is ruled out at 99.87% C.L.. Further 
we see from Tab. 2 that now (2+2) is only slightly better than (3+0) for 
LSND global. Using LSND DAR data (2+2) is somewhat less disfavoured 
with respect to (3+1). The reason for this is that for LSND DAR (3+1) gets 
much worse compared to LSND global, whereas the fit of (2+2) deteriorates 
only slightly (see Tab. 1). Hence, the difference becomes smaller. Still, (3+1) 
is clearly preferred: for 4 d.o.f. a Ax' = 11.5 implies that (2+2) is ruled out 
at 97.9% C.L. with respect to (3+1). Of course, for LSND DAR (3+0) is very 
strongly disfavoured with a Ax' = 28.5 because of the high significance of the 
oscillation signal in this data set, which cannot be accounted for in the (3+0) 
case. 



9 Conclusions 

Prompted by the recent SNO neutral current result and recent atmospheric 
data, we have re-analyzed the global status of all current neutrino oscilla- 
tion data in the framework of four-neutrino schemes, in a way similar to that 
of Ref. Our present update includes the global solar neutrino data, the 
most recent sample of atmospheric neutrino experiments, and a detailed treat- 
ment of experiments reporting no evidence (NEV) for oscillations (KARMEN, 
Bugey, CDHS as well as information from CHOOZ and Palo Verde). Besides 
the analysis of the global LSND data, we also perform an analysis based mainly 
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on the decay-at-rest (DAR) data sample, which provides a higher sensitivity to 
the oscillation signal in LSND. We have identified the parameter consistency 
(PC) test and the parameter goodness of fit (PG) as useful statistical methods 
to evaluate the compatibility of different data sets in a global analysis. 

We have found that in (2+2) mass schemes recent solar and atmospheric data 
are compatible only at the 3.5a according to the PC test, and at 4.8cr according 
to the PG method. In (3+1) mass schemes the disagreement of LSND with 
the rest of the oscillation data has been evaluated for both LSND analyses 
as a function of Aml^^-^. For the global LSND data PC is achieved only at 
2a and the PG is below 1% for all values of Aml^^j,. Using the LSND DAR 
sample the disagreement is always stronger: we find a PC only at 99% C.L. 
and for most values of Am^gj^^ the PG is worse than Aa. 

We have evaluated the g.o.f. of a global fit in terms of the PG by dividing 
the data in SOL, ATM, LSND and NEV samples. In the best case we find 
a PG of only 0.56%. This value occurs for the (3+1) mass scheme and the 
global LSND data. Using the LSND DAR analysis we get a PG of 7.6 x 10"^ 
for (3+1). For (2+2) oscillation schemes the situation is worse: we find a 
very bad PG of 1.6 x 10"^ for LSND global and 3.5 x 10"^ for LSND DAR. 
Concerning the relative status of the hypotheses (3+1), (2+2) and the three- 
active neutrino case (3+0) we find that for the LSND global data (2+2) and 
also (3+0) are strongly disfavoured with respect to (3+1) with a Ax^ = 17.8 
and Ax^ = 20.0, respectively. For LSND DAR we find for (2+2) a Ax^ = 11.5 
relative to (3+1), and the high significance of the oscillation signal condensed 
in the DAR sample leads to the huge value of Ax^ = 28.5 of (3+0) relative to 
(3+1). 

The exclusion of four-neutrino oscillation schemes of the (2+2)-type is based 
on the improved sensitivity of solar and atmospheric neutrino experiments 
to oscillations into a sterile neutrino, thanks to recent experimental data. 
This is a very robust result, independent of whether the LSND experiment is 
confirmed or disproved. The exclusion of (3+1) schemes depends somehow on 
the used LSND data sample^. Furthermore, it heavily relies on the results of 
negative SBL experiments, especially on the Bugey and CDHS disappearance 
experiments. Therefore, if LSND should be confirmed by MiniBooNE, it will be 
crucial to improve the experimental data on SBL 'z/g and/or ^ly^ disappearance. 
If the signal in LSND should indeed stem from oscillations in a (3+1) mass 
scheme the SBL and VJ, disappearance amphtude must be just on the edge 
of the sensitivity of current experiments. 

^ We want to stress that, although (3+1) models themselves are not ruled out, 
they do not offer an acceptable framework for a combined description of current 
oscillation data, if LSND is included. In contrast, any model of the (2+2)-type is 
not viable if the gap separating the two pairs of neutrino states is large compared 
to Ami^M- 
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In summary, the interpretation of the global neutrino oscillation data includ- 
ing LSND in terms of four-neutrino mass schemes - in either (3-1-1) or (2-1-2) 
realizations - is strongly disfavoured. In the best case we obtain a parameter 
g.o.f. of only 0.56%. The analysis we have presented brings the LSND anomaly 
to a theoretically more puzzling status indeed. We want to note that also in- 
troducing more sterile neutrinos participating in the oscillations is unlikely 
to substantially improve the situation [Q. If LSND should be confirmed by 
the results of the MiniBooNE experiment the situation will become even more 
puzzling. Examples of more far-fetched alternative explanations are the possi- 
bility of lepton number violating muon decay or large CPT violation in the 
neutrino sector Such scenarios will be crucially tested at the upcoming 



experiments MiniBooNE EO] and KamLAND [|| . 



A Statistical compatibility of different data sets 



In the global analysis of oscillation data in the framework of four-neutrino 
schemes we are faced with the problem that for a given mass scheme some 
data are in conflict with each other. In this appendix we present two statistical 
methods to evaluate quantitatively whether such data sets are consistent with 
each other. 

Suppose a theoretical model characterized by P parameters x = (xi, . . . ,xp) 
and K statistically independent experiments (or data sets) described by the 
theory. Each of these experiments may consist of many data points; experiment 
k observes A^^ data points. Furthermore, experiment k may in general depend 
only on Pfc < P parameters , which is a sub-set of all parameters x of the 
full theory, leading to the x^-function 

xl{x') = {xlUn + ^xl{x'). (A.l) 

The three terms are distributed with A^, Nk — Pk, Pk degrees of freedom 
(d.o.f.). The global is given by 

XL(^) = E Xl{x') = xLn + ^X\x) . (A.2) 
fc=i 

Usually the goodness-of-fit (g.o.f.) is evaluated by considering the value of 
Xmin' which is distributed with N — P d.o.f., where A^ = X)fc -^fc- 

Let us assume, that all the experiments considered separately lead to a good 
fit: (Xfc)mm ^ Nk — Pk- However, there may be some disagreement between 
different experiments. If the number of data points is very large the disagree- 
ment may be completely washed out and the usual g.o.f. would be very good. 
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i.e. Xmin ^ N — P despite the disagreement .Q Here we propose two methods 
which are more useful in global analyses with many data points, like in the 
four- neutrino case. 



A.l Parameter consistency test (PC) 



Let us consider the case of two data sets {K = 2), which have P < P parame- 
ters X in common. Then at a given CL. the allowed regions for the parameters 
X according to each data set are determined by two functions Axi{x) and 
Ax2{x), where we minimize with respect to the other parameters. Now we cal- 
culate Xpc; defined as the minimal Ax^ value where there is still some overlap 
between the allowed parameter regions of the two data sets. This Xpc ^as to 
be evaluated for P d.o.f. to obtain the CL. at which there exists still some 
values of the parameters such that both data sets are consistent. 



We call this criterion parameter consistency test. It can be applied if the global 
data can be divided into two potentially conflicting sets. Similar methods have 
been used in Refs. 



36], ^ to evaluate the consistency of LSND and NEV 
experiments in (3+1) schemes and in Ref. [0 to check the compatibility of 
LSND and KARMEN. 



A. 2 Parameter goodness of fit (PG) 



Let us write the global x^ Eq. (A. 2) in the following way: 



XIM) = E(xDmin + E ^Xlix") . (A.3) 
k=l k=l 

To evaluate the g.o.f. we now use only the second term: 

fix) ^ f: Axlix') = xLn + Ax'(x) • (A.4) 
fe=i 

The first term in Eq. (A.3) describes the g.o.f. of all the experiments taken 
alone, whereas xLin Eq. (A.4) corresponds to the contribution to the total 
minimum Xmm which stems from the combination of the data sets. The three 
terms in Eq. (A.4) are distributed with J2k Pk, J2k Pk — Pi P d.o.f., respectively. 



' For the global four-neutrino analysis this effect was discussed explicitly in 
Ref. [|8|. 
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Hence, to check the g.o.f. we just have to evaluate 

XPG = xLn for J^Pk-P d.o.f.. (A.5) 

k=l 

We call this criterion parameter goodness of fit and in general it is more 
stringent than the parameter consistency test (for an illustration see Fig. 6). 

If Mp < K experiments depend on the parameter Xp we have the obvious rela- 
tion Ylik=i Pk — Y^p=i Mp, and the relevant number of d.o.f. given in Eq. (A.5) 
can be written as J2p=i Mp — P = J2p=iiMp — l). From this relation we see that 
parameters relevant for only one experiment {Mp = 1) do not contribute be- 
cause they are minimized and contribute only to the individual (xDmin which 
are not considered in Xpg- What is tested is only the parameter dependence 
of the data sets relevant for the combination. 

A simple example: Let us consider the case of only one parameter y, which is 
determined by 2 experiments. Hence, we have K — 2 and P — Pi — P2 — 1 
and XpG — n^i^ [^Xi(y) + ^xK?/)]' which has to be evaluated for 1 d.o.f.. The 
following statement is easy to prove: Let yi and y2 be the best fit values of the 
two experiments (Axliyk) = 0). Then the C.L. at which yi — y2 is consistent 
with zero is identical to the parameter goodness of fit. 
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